Exact Solution of a Charge- Asymmetric 
Two-Dimensional Coulomb Gas 

L. Samaj 1 
February 1, 2008 

Abstract 

The model under consideration is an asymmetric two-dimensional Coulomb gas 
of positively (qi = +1) and negatively (52 = —1/2) charged pointlike particles, 
interacting via a logarithmic potential. This continuous system is stable against 
collapse of positive-negative pairs of charges for the dimensionless coupling constant 
(inverse temperature) (3 < 4. The mapping of the Coulomb gas is made onto the 
complex Bullough-Dodd model, and recent results about that integrable 2D field 
theory are used. The mapping provides the full thermodynamics (the free energy, 
the internal energy, the specific heat) and the large-distance asymptotics of the 
particle correlation functions, in the whole stability regime of the plasma. The 
results are checked by a small-/? expansion and close to the collapse (3 = 4 point. 
The comparison is made with the exactly solvable symmetric version of the model 
(gi = +1,^2 = — 1); and some fundamental changes in statistics caused by the 
charge asymmetry are pointed out. 
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1 Introduction 



A two-dimensional (2D) Coulomb plasma is the continuous system of charged particles, 
plus perhaps a uniformly charged neutralizing background, confined to a plane and in- 
teracting via the logarithmic Coulomb potential. In this paper, the classical equilibrium 
statistical mechanics at the (dimensionless) inverse temperature (3 is studied. We will 
restrict ourselves to the simple point-particle Coulomb systems: if there are at least two 
species of charged particles with opposite signs, the Coulomb system is stable against 
the charge collapse for small enough (3 < /3 co i (in 3D, the collapse point /3 co i — > and 
the stability requires quantum mechanics [1]). Going beyond the collapse point needs a 
short- distance (e.g. hard core) regularization of the Coulomb interaction. The famous 
Kosterlitz-Thouless (KT) transition [2] of infinite order from a high-temperature conduct- 
ing phase to a low-temperature insulating phase takes place at some Pkt > Pco\- Like 
in other dimensions, the 2D Coulomb systems admit the Debye-Hiickel treatment in the 
high-temperature limit (3 — > 0, as was rigorously proven in Ref. [3]. As concerns the bulk 
thermodynamics, based on the scale invariance of the logarithmic potential, the density 
derivatives of the free energy, like the pressure, can be obtained exactly in the whole 
stability range of temperatures [4]- [6]. On the other hand, the temperature derivatives 
of the free energy, like the internal energy or the specific heat, are nontrivial statistical 
quantities which cannot be obtained by simple means and they are known only in special 
cases (see below). 

The most studied versions of the Coulomb plasma are the one-component plasma 
(OCP) and the symmetric two-component plasma (TCP), or Coulomb gas. 

The 2D OCP of equally, say unit, charged pointlike particles in a uniform neutralizing 
background is formally related to the fractional quantum Hall effect [7, 8]. No collapse 
occurs. There are indications that, around (3 ~ 142, the fluid system undergoes a phase 
transition to a 2D Wigner crystal [9]. In a more recent paper [10], the existence of this 
transition has been put in doubt. The model is exactly solvable at (3 = 2 by mapping onto 
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free fermions, in the bulk [11] as well as in some inhomogeneous situations (see review 
[12])- 

The symmetric 2D TCP consists of oppositely ±1 charged particles, no background is 
present. The collapse point f3 co \ = 2 coincides with the exactly solvable free-fermion point 
of the equivalent Thirring model [13]- [15]. At this collapse point, for a fixed fugacity, 
while the free energy diverges, truncated Ursell functions are nonzero and finite. Quite 
recently [16], the complete bulk thermodynamics of the symmetric 2D TCP was derived 
exactly in the whole stability region j3 < 2. The mapping onto a bulk 2D sine-Gordon 
theory with a conformal normalization of the cos-field was made, and recent results about 
that field theory were applied. Subsequently, the surface tension of the same model in 
contact with an ideal-conductor [17] and with an ideal-dielectric [18] rectilinear walls was 
obtained via a mapping onto integrable 2D boundary sine-Gordon theories with Dirichlet 
and Neumann boundary conditions, respectively. The large-distance behavior of the bulk 
charge [19] and density [20] correlation functions was derived exactly by exploring the 
form-factor theory of the equivalent sine-Gordon model. 

Hansen and Viot [21] have introduced a general asymmetric 2D Coulomb gas which 
consists of two species of pointlike particles with positive and negative charges of arbitrary 
strengths q a (a = 1,2). For [gi/^l being an integer, the model describes a plasma of 
electrons and ions of integer valence. Without any loss of generality, we choose 

gi = +1 and q 2 = -1/Q (1.1) 

where Q — 1, 2, . . . is a positive integer. The model interpolates between the symmetric 
TCP (Q = 1) and, after subtracting the kinetic energy of 2-species, the OCP obtained 
as the extreme asymmetry case Q — > oo. The Boltzmann factor of a pair of positive 
and negative charges, r _/3//Q , is integrable at short distance in 2D space for (3 < 2Q, so 
the collapse point is /3 co i = 2Q [21]. In the case of a vanishing but nonzero hard core 
around particles, the KT phase transition was conjectured to take place at /3kt = 
[22]. Highly asymmetric Coulomb mixtures in the strong-coupling regime have attracted 
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much attention in the last years [23] -[26] due to the phenomenon of overcharging (charge 
inversion), i.e. the situation when the number of counterions in the vicinity of a macroion 
is so high that the macro-charge is overcompensated. 

In this work, we report the exact solution of the Q = 2 asymmetric 2D Coulomb 
gas, with the lowest degree of charge asymmetry q 1 — 1 and q 2 = —1/2, in the whole 
stability interval of pointlike particles (3 < 4. The exact solution involves the complete 
bulk thermodynamics (the free energy, the internal energy, the specific heat) and the 
large-distance behavior of the particle correlation functions. It is obtained by mapping 
the underlying Coulomb gas onto a member of 2D Toda field theories [27], namely the 
complex Bullough-Dodd (cBD) model [28]-[30] with a conformal normalization of the ex- 
ponential field. Recent results about that integrable field theory, derived by using the 
Thermodynamic Bethe Ansatz (thermodynamics) and the form-factor method (correla- 
tion functions), are applied. Since the calculations in the cBD model are based on special 
analyticity assumptions, which are not yet rigorously proven, we check the results for the 
plasma by a small-/? expansion, using a renormalized Mayer expansion in density [31, 32], 
and close to the collapse (5 — 4 point, using an electroneutrality sum rule [33] combined 
with an independent-pair conjecture made by Hauge and Hemmer [6]. The comparison is 
made with the symmetric version of the model, and the fundamental changes in statistics 
due to the charge asymmetry are pointed out. 

The paper is organized as follows. In Section 2, all relevant aspects of the mapping 
between the asymmetric 1/ — \ 2D Coulomb gas and the cBD model are presented. The 
complete thermodynamics of the Coulomb system is derived in Section 3. The asymptotic 
large-distance behavior of the particle correlation functions is obtained within the form- 
factor method in Section 4. Section 5 is a brief recapitulation with some concluding 
remarks. The results are checked by the small-/? expansions in Appendix A and close to 
the collapse (3 — 4 point in Appendix B. 
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2 Field-theoretical representation of the asymmetric 
2D Coulomb gas 

We consider the asymmetric TCP made up of two species of pointlike particles with 
charges q a (a = 1, 2) given by (1.1). The particles are confined to an infinite 2D space of 
points r e R 2 , and interact with each other by the pair Coulomb potential. The Coulomb 
potential v at spatial position r, induced by a unit charge at the origin, is given by the 
2D Poisson equation 

Av(r) = -27r5(r) (2.1) 

as follows 

V (r) = -ln(|r|/r ) (2.2) 
The length constant r , which fixes the zero point of the Coulomb potential, will be set 
for simplicity to unity. The interaction energy E of particles {j} reads 

£ = E^^(|r,-r fe |) (2.3) 

j<k 

Introducing the microscopic density of particles of species a, n a {r) = J2j^a,aj^( r — 
the microscopic densities of the total particle number and of the charge are 

n(r) = £ n a (r), p(r) = £ g CT n CT (r) (2.4) 

respectively. The energy (2.3) can be thus written as 

E=Ud 2 rf dVp(rM|r - r'|)p(r') - i„(0) £ <. (2.5) 
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where v(0) is the (divergent) self-energy. 

We will work in the grand canonical ensemble, with position-dependent fugacities z a (r) 
of species a — 1,2. The grand partition function S at inverse temperature f3, considered 
as a functional of the species fugacities, is defined by 

S[^,^ 2 ] = E E ¥T¥T/ II [ d \^i( u i) 

Ar 1=0 JV 2 =0 7V 1 ! 7V 2 ! J j=i 

N 2 

x II [^-^(Vj)] ex P [-P e n u n 2 {u,v)\ (2.6) 

J'=l 



where En 1: n 2 (u,v) denotes the Coulomb interaction energy (2.3), resp. (2.5), of Ni parti- 
cles of type 1 at vector positions {u J }^ 1 and N 2 particles of type 2 at positions {vj}^l v 
The statistical quantity S can be expressed in terms of a 2D Euclidean theory with the 
aid of the standard procedure (see, e.g., ref. [35]). Using the representation (2.5) for E in 
exp(— f3E) and assuming that — A/(27r) is the inverse operator of v(r) [see relation (2.1)], 
one applies the Hubbard-Stratonovich transformation 



exp 



^Jd 2 rjd 2 r'p(r)v(\r-r'\)p(r') 



/ Vcj) exp 


'/ d 2 r 




/P0exp (jdV^A^ 





(2.7) 



where 4>(r) is a real scalar field and / V(j) denotes the functional integration over this field. 
The term 0A0 can be turned into — (V0) 2 by using integration by parts. The summation 
over Ni and N 2 in (2.6) then implies 

/£>(/> exp (S[z 1 ,z 2 \) 



^ Zl,Z2] /2>^exp(-5[0,0]) 
where the action of the 2D Euclidean field theory takes the form 

1 



S[z 1 ,z 2 ] 
b 2 



d 2 r 

(3/A 



16tt 



zi r e 



z 2 (r)e 



■i(b/Q)4> 



(2.8) 

(2.9a) 
(2.9b) 



Here, the fugacities z a (a = 1, 2) are renormalized by the self-energy terms exp[/3v(0)q 2 ./2]. 
For the homogeneous system with uniform species fugacities, z a (r) = z a , the action (2.9) 
simplifies to 



S(z u z 2 ) = jd 2 r 



16tt 



(V0) 5 



Z\Q 



z 2 e 



-i(b/Q)4> 



(2.10) 



Note that this action is complex, except in the symmetric case Q — 1 with z\ — z 2 . 

The field representation of the multi-particle densities can be obtained from the func- 
tional generator S, defined by (2.8) and (2.9), in a straightforward way. At the one-particle 
level, the density of particles of type a(— 1, 2) is given by 



n n 



(Mr)> 
1 5E 



ESz a (r) 
6 



uniform 



(2.11) 



Consequently, 

n a = z a (e ib ^) (2.12) 

where (• • •) denotes the averaging over the action (2.10). Note that under the shift of 
the field variable — > <p + <p in (2.10) (which has no effect on S as the whole), the 
species fugacities change as Z\ — > Ziexp(i60 o ) an d z 2 — > -22 exp[— i{b/Q)4>o\, and therefore 
S depends only on the combination Zi(z 2 ) Q . According to the definition (2.11), this 
property of S implies the neutrality condition n\ — (1/Q)n 2 = 0. At the two-particle 
level, one introduces the two-body densities 



?w(r, r') = (n a (r)n a/ (r')) - nj aa/ 5(r - r') 
1 S 2 E 



uniform 



(2.13) 



E5z a (r)5z al (r') 
so that 

n aa ,(\r - r'|) = z.vie^e^'^) (2.14) 

It is useful to introduce also the pair distribution functions 

, rw(r,r') 
g cc , (r , r ) = (2.15) 

and the (truncated) pair correlation functions 

K al {v,v')=g aal {v,v')-l (2.16) 

In statistical mechanics, the short- distance behavior of the pair distribution function 
is dominated by the Boltzmann factor of the pair Coulomb potential [21, 34], 

W(r,r') ~C^|r-r'|^' as |r - r'| -> (2.17) 

(provided /3 is small enough as explained below). The prefactors C aa i are expressible as 
follows 

C aa , = exp [p + »~ - »l +qa )] (2.18) 
where /i^ is the excess chemical potential of species a defined by 

='«(!) < 2 - i9 > 
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and with q arbitrary represents an extended definition of the excess chemical potential 
for a "guest" particle of charge q put into the considered Coulomb gas (1.1). To get the 
stability region of /i^ x , one has to consider the interaction Boltzmann factor of the g-charge 
with an opposite charge from the plasma. When q > (q < 0), the Boltzmann factor with 
the opposite —1/Q (+1) plasma charge at distance r, r _/3<? / Q (r l3q ), is integrable at small 
r if and only if (3 < 2Q/q {(3 < -2/q). ^ is therefore finite if (3 < (Q + l)/\q\ + (Q-l)/q 
and goes to — oo outside of this stability region. According to the formula (2.18) this 
means that the prefactor C i2 remains finite in the whole stability region j3 < 2Q, while 
the prefactors C\\ and C 22 are finite only if (3 < Q. The divergence of Cn and C 2 2 in 
the middle of the stability region, at the point (3 = Q, is accompanied by a change in 
the short-distance behavior (2.17) of g u and g 22 . The general analysis in Ref. [21] shows 
that, in the short- distance limit r — > 0, 

9n(r) oc r' 3 , (3 < Q 

r 2m+f3[l-m(4Q+l-m)/(2Q 2 )] 2 Q 2 n 2 Q 2 2Q x 

2Q-m+l 1 2Q-m v ' ; 
where m is an integer from the interval 1 < m < Q, and 

g 22 (r) oc r^ 2 , (3 < Q 

oc r 2-/3(2Q-i)/Q 2 ) Q<(3<2Q (2.20b) 

Close to the collapse point (3 — * 2Q, one finds 

9n (r) oc r®- 1 (2.21a) 
522 (r) oc r i2 ' Q ^ 2 (2.21b) 

The collapse relation (2.21b) tells us that, for Q > 1, g 22 (r) does not vanish as r — > 
(as is intuitively expected since equally charged particles repel each other), but diverges. 
This is related to a paradoxical clustering of counterions in the asymmetric plasma. 
We now intend to derive the short-distance behavior (2.17), say for g u , in terms of 
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the 2D field theory (2.10). According to the definitions (2.6), (2.13) and (2.15), we have 
<?n(r,r) = \ v -T\ p [ — \ - £ E N \ N \ / II [ d u M^j)\ U d ^(v,-) 

1 1—1 Ni=2 N2=0 Jv l- iy 2- J j=i J= i 

xexpj-^^K^)-/? | dV'Kr-r")+f(r'- r")]p(r")} (2.22) 

where p denotes the microscopic charge density of Aq — 2 particles of type 1 and N 2 
particles of type 2. One shifts Aq — > iVi — 2 in the summation over Aq. When r' — > r, the 
exponential in (2.22) can be written in the form presented on the lhs of Eq. (2.7) with 
the substitution p(r") — > p(r") + 25(r" — r). Using the Hubbard-Stratonovich identity 
(2.7), the previously generated single term i^J (3/A <pp now involves two terms: (3/A <pp + 
i2y / ' (3 1 '40(r) . One integrates over particle coordinates, except the fixed r-coordinate, and 
ends up with 

Sudr-r'l) ~ \v-r'f ^^( e i2 ^) as |r - r'| -> (2.23a) 

Analogously, one finds that 

<7i 2 (|r-r'|) ~ Ir-r'j-^^^^eKi-VQ)^) as | r _ r '|^ (2.23b) 

<?22(|r-r'|) ~ \ Y - r f/Q 2 (^l\\ e -^/Q)4>) as | r _ r '|^ (2.23c) 

\ri2J 

Comparing these relations with Eqs. (2.17) - (2.19), the excess chemical potential of a 
particle with charge q embedded into the considered plasma is related to the one-point 
expectation of the exponential field as follows 

cxp (-(3pf) = (e 1 ^) (2.24) 

Note that relations (2.12) are special cases of the general formula (2.24). According to 
the analysis after Eq. (2.19) and with respect to the relationship (2.9b), the one-point 
expectation (e iqb ^) is finite if 4b 2 < (Q + l)/\q\ + (Q — 1)/?, and goes to +00 otherwise. 
Comparing relations (2.23) with Eqs. (2.14) and (2.15), the short-distance behavior of 
two-point expectations of the exponential field is given by 

^(r) e i g W)} „ | r _ r '|4& 2 ?? '( e i( g+9 ')^} ag | r _ r '| ^ (2.25) 
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In quantum field theory, the model (2.10) can be regarded as the Gaussian confor- 
mal field theory perturbed by the operator —z\ exp(i&0) — z 2 exp[— i(b/Q)<f>]. The species 
fugacities z\ and z 2 are renormalized by the (divergent) self-energy factors. To give the 
parameters z\ and z 2 a precise meaning, one has to fix the normalization of the adjoint 
exponential fields. The normalization which corresponds to the short- distance limit of 
the two-point function (2.25) is known as conformal. Under the conformal normalization, 
the diverging self-energy factor disappears from statistical relationships. This makes the 
bridge between the underlying asymmetric Coulomb plasma and the corresponding field 
theory with the action (2.10), treated within the Conformal Perturbation theory. 

Note that q and q' should be sufficiently small to ensure that (2.25) is the leading short- 
distance asymptotics. From the point of view of the above equivalence, the crucial is the 
stability of the two-point expectation asymptotics (2.25) for two opposite charges q — 1 
and q' = -1/Q. In that case, the stability region of (e i ^ +9 ') 6 ^) is (3 < /3 stab = 2Q 2 /(Q- 1); 
for the symmetric Coulomb gas /3 sta b - ► °o, for Q = 2 /3 stab incidently coincides with 
Pkt — 8 an d for any finite Q it holds /3 sta b > f3 co \. 



3 Thermodynamics of an asymmetric Coulomb gas 



There is a large class of integrable 2D field theories, known as the afnne Toda theories, 
which are based on the Dynkin-diagram classification of simple Lie groups (for a nice 
review, see Ref. [27]). Let Q be a simple Lie algebra of rank r, {e 1; • • •, e r } a set of 
simple r-dimensional roots of the corresponding Dynkin diagram and — e = J2l=i n i e i 
(the coefficients {n,j} are called Kac labels) the maximal root. The affine Toda theory 
built on Q is defined by the action 



S 



J 



d 2 r 



^ (d,<f>) 2 + £ Zie bei -+ + z r+1 e be -<t> 



i=i 



(3-1) 
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where the field 4> = (0i, ■ • • ,cf> r ) consists of r real scalar components and b is the real 
coupling constant. The fields in Eq. (3.1) are normalized so that at {zi = O}^ 1 

(0 a (r)0 fe (r / )) = -45 a5 ln|r-r'| (3.2) 

For the one-component (r = 1) case of interest, there exist two integrable Toda theories. 
When Q = Lie group with e± — 1 and e = — 1, one obtains the sinh-Gordon (b real) 
or sine-Gordon {b pure imaginary) models. The sine-Gordon model is identified with the 
action (2.10) with Q — 1 and describes the thermodynamics of the symmetric 2D TCP 
(see Refs. [16]-[19]). When Q = Lie group with e\ — 1 and e = —1/2, the action 
(3.1) takes the form 



>BD 



f d 2 r [— (V0) 2 + ztf* + z 2 e~ b ^ 2 



(3.3) 



This theory is known as the Bullough-Dodd (BD) model [28]. Its complex version, ob- 
tained via the substitutions £» — > i£» and z± j2 — > —^1,2, 



^cBD — 



[ d 2 r \-L (V0) 2 - z x <P* - z 2 e~^ 



(3.4) 



is referred to as the Zhiber-Mikhailov-Shabat model [29, 30], or simply the complex 
Bullough-Dodd (cBD) model. Comparing (3.4) with (2.10) one observes that the cBD 
model with the short- distance normalization (2.25) is the 2D field realization of the Q = 2 
asymmetric TCP, i.e. the system of +1 and —1/2 charged particles. According to (2.9), 
the coupling constant b is related to the inverse temperature by b 2 = (3/ 4 and the param- 
eters Z(j (a = 1,2) represent the renormalized species fugacities. 

The particle spectrum of the BD model (3.3) consists of a single neutral particle of 
mass to. The two-particle ^-matrix was described in Ref. [36]. The (dimensionless) 
specific grand potential 

— uj = lim ^rlnS (3.5) 

V^oo V 

where S is given by (2.8) and V is the volume, was obtained in terms of the particle mass 
to in Ref. [37] following the Thermodynamic Bethe Ansatz technique [38, 39]. In the 
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same Ref. [37], the relation between the particle mass m and the model parameters zip 
was established under the conformal normalization, and an explicit formula for the mean 
value of the exponential field (e a< ^) was suggested by exploring a reflection relationship 
between the BD model and the 2D Liouville theory [40, 41]. The derivation of the results 
was based on special analyticity assumptions, so their verification by various checks is 
needed. 

The spectrum of the cBD model (3.4) exhibits an extremely complicated hierarchy of 
particles [42]. The fundamental particle is a three-component kink. The kinks generate 
one-component bound states (breathers) and higher kinks, these higher kinks generate 
new breathers and new kinks, etc. The important simplifying fact is that, in the whole 
stability interval of interest < b 2 < 1 (0 < (3 < 4), the lightest particle, the 1-breather, 
corresponds to the analytic continuation of the only particle in the spectrum of the BD 
model. Since the lightest particle dominates in the thermodynamic limit V — > oo, we 
can apply all results of Ref. [37] presented in the above paragraph, with the substitutions 
£» — > i£» and z ij2 — > —^1,2, also to the cBD model. In particular, the specific grand potential 
(3.5) takes the form 

2 

m 

—uj = = (3.6a) 

16v^sin(7r£/3)sin(7r(l + 0/3) 



b 



,2 



2-b 2 

and the mass of the lightest 1-breather reads 



(3.6b) 



2 v / 3r(l/3) 

m = 



r(i-e/3)r((i + o/3) 



Zl nT(l - 6 2 )] {1+ « /6 [ 2z 2 7rr(l-&V4) l (] " ; ' 

r(6 2 /4) 



(3.7) 



T(b 2 ) 

Note that at the collapse point b 2 = 1 (f3 = 4) the particle mass m — > 00. The expected 
divergence of — uj (3.6) at b 2 = 1 is therefore caused by m 2 , while the prefactor to m 2 
remains finite. Such behavior is opposite to that observed in the symmetric 2D TCP [16] 
at the corresponding collapse point where the mass of the lightest sine-Gordon breather, 
mi, is finite and — uj diverges due to the prefactor to m\. The expectation value of the 
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exponential field is given by (see also [43]) 

~z 2 2- b2 ' 2 T(l + 6 2 )r(l - 6 2 /4) 



(e ia0 > 



zi r(i -6 2 )r(i + 6 2 /4) 



mT(l-e/3)r((l + 0/3) 



2 2 / 3 v / 3r(l/3) 



2a 2 -ab 



dt 

x< ' xl) Wo T 



sinh((2 -6 2 )t)^(t,a) 



sinh(3(2 - 6 2 )t) sinh(2t) sinh(6 2 t) 



2a z e 



2 -It 



(3.8a) 



where 



V(t,a) = -sinh(2afa) 



sinh((4 - b 2 - 2ab)t) - sinh((2 - 26 2 + 2ab)t) 
+ sinh((2 - b 2 - 2ab)t) - sinh((2 - b 2 + 2a6)f) 
-sinh((2 + 6 2 -2ab)t) 



(3.8b) 



The integral in (3.8a) is well defined if 



-1 < Re(a) < 1 



(3.9) 



Taking into account (2.12), the charge- neutrality condition in the considered Coulomb 
gas, ni = TI2/2, results in the equality 



Zl (e#)^( e -W^) 



Using the integral formula for the logarithm of the Gamma function [44] 



(3.10) 



\nT(z 



. r°° dt 

)= i r 



~t 



(z-l) 



"(*-!)* _ 1 



Re{z) > 



(3.11) 



1 -e"* 

it can be readily verified that the suggested formula (3.8) is consistent with the neutrality 
condition (3.10). 

We are now ready to derive the basic density-fugacity relationship for the asymmetric 
1/ — I Coulomb gas. The total particle number density is 

d(-u) d(-cu) 



n — zi- 



+ z 2 - 



= (i + 0(-") 



(3.12) 



dzi dz 2 

where the auxiliary parameter £, introduced in (3.6b), is expressible in the inverse tem- 
perature (3 = Ab 2 as follows 



g 
8-/3 



(3.13) 
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After some algebra, Eqs. (3.6) and (3.7) give the explicit density-fugacity relationship 



n 



1-/3/8 



(^2 2 ) 1/3 



Vs r 2 (i/3) r(£/3)r((2-o/3) 



x 



1 (1 - /3/8)tt 2 


r(i 


-C/3)r((i + 0/3) 




1/3 


"27rr(l -/3/16)" 


2/3 


[ r(/3/4) . 




[ r(/3/i6) 





1-/3/8 



(3.14) 



The length constant r in (2.2) was set to unity. This can be shown to imply that the 
fugacity product Z\z\ in neutral configurations of (2.6) has dimension [length] -st 1 -/ 3 / 8 )^ 
and so Eq. (3.14) is dimensionally correct. The small-/? expansion of the rhs of (3.14) 
reads 



n 



1-/3/8 



exp 



3^(2/3) - 2tt 2 



P 2 



1728 



+ 0(/3 3 ) (3.15) 



where C is the Euler number and ip(x) = d[lnT(a;)]/da; is the psi function. The series 
representation of the first derivative of tp reads 

1 



(3.16) 



The series expansion (3.15) is checked up to the indicated order in Appendix A.l by using 
a bond-renormalized Mayer expansion in density [31, 32]. For a fixed fugacity product 
ziz%, relation (3.14) implies the expected collapse singularity of the density n: 



n 



(^1) 2/3 



2tt 



r(3/4) 



14/3 



1 



(l-/3/4) 2 /3 



as (3 



(3.17) 



r(i/4)_ 

This singularity is reproduced indirectly in Appendix B.l by applying a (slightly modified) 
perfect-screening sum rule, which is another important check of the basic result (3.14). 

To obtain the complete thermodynamics of the asymmetric 1/ — \ Coulomb gas, we 
pass from the grandcanonical to the canonical ensemble via the Legendre transformation 



F{T; Ni, N 2 ) = Q + ^ + fi 2 N 2 



(3.18) 



where 



-u((3,n) 



k B Tu((3,n)V 
(-> ^ 



8 



n 



(3.19a) 
(3.19b) 
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Ni = N/3, N2 = 2N/3 with N = nV being the total particle number, and 



Ha((3,n) = k B T \nz a ((3,n), (a = 1,2) 



(3.20) 



is the chemical potential of species a. The dimensionless specific free energy, / = 
Fj (NkuT) , is thence given by 



(3.21) 



With the aid of the density-fugacity relationship (3.14), one has explicitly 



r 2 (i/3) r(e/3)r((2-p/3) 
(i-/3/8)r(i-e/3)r((i + o/3) 



r(i-/9/4) 
r(/3/4) 



3 



2r(l -/3/16) 



r(/3/i6) 



(3.22) 



The (excess) internal energy per particle, u ex = (E)/N, and the (excess) specific heat at 
constant volume per particle, Cy = Cy/N, are determined by elementary thermodynam- 
ics as follows 



(3.23a) 
(3.23b) 



The specific heat is independent of the particle number density n, which is a specificity 
of the 2D pointlike Coulomb gases. The expansion of Cy near the collapse point (3 = A 
results into the Laurent series 



„ex 
Cy 



+ 0(1), /3-4" 



(3.24) 



k B 3(1 -/3/4) 2 3(1-/5/4) 

The leading two terms in (3.24) are reproduced indirectly in Appendix B.2 by applying 
an independent-pair conjecture of Hauge and Hemmer [6]. 
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4 Large-distance behavior of particle correlations 

In a 2D integrable field theory with particle spectrum {e}, correlation functions of local 
operators O a (a is a free parameter) are expressible as infinite convergent series over multi- 
particle intermediate states (see e.g. Ref. [45]). For the truncated two-point correlation 
functions 

(O a (r)O a ,(r')) T = {O a {v)O a ,{v')) - {O a ){O a ,) (4.1) 

the series reads 

(O a (r)O a ,(r')h = E ) 9 '^ N F a (9 1: ...,9 N ) €l ... eN 

N=l iV! ei,...,^- 7 " 00 { Z7T ) 

x eN -^F a ,(9 N , ...,9 ± ) exp ^-|r - r'| £> e . cosh^j (4.2) 

Here, 9 e (— oo, oo) is the rapidity which parametrizes the energy E and the momentum 
p of a particle e of mass m e as follows 

_E = m e cosh^, j9 = m e sinh^ (4.3) 

and the normalization constants in the form factors {F a } depend on the specific form of 
the operator O a . In what follows, we will consider O a to be an exponential field: 

a(r) = exp(a0(r)) (4.4a) 

in the case of the BD model with the action (3.3), and 

O a (r) = exp (ia0(r)) (4.4b) 

in the case of the cBD model with the action (3.4). 

The form-factor representation (4.2) is particularly useful in the large-distance limit 
|r — r'| — > oo. The dominant contribution in this asymptotic limit comes from a one- 
particle intermediate state with the minimum value of the particle mass m, at the point 
of the vanishing rapidity 9 — > 0. The consequent exponential decay exp(— m\r — r'|) 
is multiplied by a slower (inverse power law) decaying function, whose particular form 

16 



depends on the one-particle form factors. For the BD model (3.3) with the only particle 
(e = 1) in the spectrum, the one-particle form factors F a {6)\ = F a and l F a i{6) = F*, are 
presented for the exponential field (4.4a) in Refs. [46, 47]. The transition to the cBD 
model (3.4) is straightforward since, as was already mentioned in Section 3, the analytic 
continuation of the single particle in the BD spectrum corresponds to the lightest particle 
of interest in the cBD model, the 1-breather with mass m given by (3.7). In particular, 
after the substitutions a — > ia and b — > ib in the one-particle form-factor formulae for the 
BD model [46, 47], the lightest-particle form factor of the exponential field (4.4b) in the 
cBD model is given by 

F a . /2-7ra A /7r / 



where 



sin(7r/3) 



1/2 



sin(27r£/3)sin(27r(l + 0/3)_ 

'_ 2 f°° dt cosh(t/6) sinh(te/3) sinh(t(l + fl/3) 1 

7o i sinhtcosh(t/2) J 1 ' J 



x exp 



According to formulae (2.12)-(2.16) adapted to the considered asymmetric 1/ — \ 
2D Coulomb gas, the pair correlation function h aa >(r,r') of species a and a' with the 
corresponding charges q a and q a i G {1, —1/2} is expressible in terms of the averages over 
the equivalent cBD model as follows 

/ ig CT 60(r) iq a ,b<t>{r')\ 

h *a' (r, r') = K . . J, . h J T (4.6) 

' / e iqabcf>\ / e iq a ,bcl>\ v > 

Using the form-factor representation (4.2) for O a (r) = exp(ia0(r)), the leading contribu- 
tion to h aa /(r,r') in the limit |r — r'| — > oo is determined by the lightest particle in the 
cBD model, the 1-breather with mass m given by (3.7) and the form factor F a given by 
(4.5). Consequently, 

F u F* l 1 foo Af) 
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Since 



/_ 



°° d6 



— mr cos 



( 7T \V2 

K (mr) ~ I J exp(— mr) 



(4.8) 



at asymptotically large r (fT is the modified Bessel function of second kind), we finally 
arrive at 

x 1/2 



h aa >(r) ~ -A^ ^— J exp(-mr) 



(4.9) 



Here, 



8^ 



7T 



exp 



00 dt cosh(t/6) sinh(^/3) sinh(t(l + f ) /3) 1 
t sinhtcosh(t/2) J 

J /2?r . \ /2tt 



" sin(2</3) sin(2 7 r(iT07^) ^ (iT* 7 *) ^ (l^ 
x cos (1 + 2£ - 4£ g<r )) cos (| (1 + 2£ - 4^)) 



(4.10) 



and the mass m (3.7) is expressible in terms of the inverse Debye length for the considered 
Coulomb gas 

K=(7r(3n) 1/2 (4.11) 



by combining Eqs. (3.6) and (3.12), 

'2V3 . /vr£ 



171 — K, 



^ 8111 ItJ 8111 U (1+0 



1/2 



(4.12) 



We recall that f = /?/(8 - 

In the high-temperature limit, the parameters \ aa > (4.10) and m (4.12) of the large- 
distance asymptotics (4.9) have the following small-/? expansions 

' 1 



A„ 



7T 



+ — -= [6(q a + q c ,) - 1] + 0(/3 3 ) (4.13a) 



24 72^ 



m 



K 



l + + 0{(3 2 ) 

48^ 



(4.13b) 



These small-/? expansions are checked up to the indicated order in Appendix A. 2 by using 
the renormalized Mayer expansion. The leading-order terms \ aa > ~ I3q a q a ' and m ~ k 
correspond to the Debye-Huckel approximation, the first corrections to this approximation 
are implied by the renormalized Meeron graph. 
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Let us introduce the charge density and particle number density combinations of the 
pair correlation functions 

h p (r,r') = ^n a q a n al q al h aa ,(j,r') (4.14a) 

a, a' 

h n (r,r') = ^n a n a ,h aa/ (r,r') (4.14b) 

a, a' 

respectively. It is evident that within the Debye-Hiickel approximation the two-particle 
correlations are determined at large distances by the charge-charge correlation function 
hp, while h n is identically equal to zero. Taking into account the first correction in X aa i 
(4.13a), the strong division between h p (r) and h n (r) disappears: both of them decay at 
large r exponentially with the same correlation length = 1/m, only the corresponding f3- 
dependent prefactors differ from one another. The prefactor is of the form X aa i = q a q a i\((3) 
in the case of the symmetric 2D TCP [19, 20], and so at any (3 in the stability regime 
the large-distance asymptotics of the two-particle correlations are determined exclusively 
by hp (h n is related to the heavier 2-breather and therefore goes to zero faster then h p ). 
The asymmetry in the particle charges thus causes a fundamental change in the relative 
large-distance behavior of the charge and density correlation functions. 

Since k — > oo at the collapse point (5 — 4, the particle mass m (4.12) diverges, and 
h aa i given by (4.9) reduces trivially to zero. On the other hand, the mass of the lightest 
particle is finite (for a fixed z) at the collapse point for the symmetric 2D TCP [19], 
and the corresponding Ursell functions U aa '(r, r') = n a n a 'h aa '(r, r') have a nontrivial 
large-distance dependence. Also from this point of view the asymmetry in the particle 
charges has a relevant influence on the large-distance characteristics of particle correlation 
functions. 

5 Conclusion 

In this paper, we have solved exactly the 2D Coulomb gas of pointlike charged particles, 
with the charge asymmetry q± — 1 and q 2 = —1/2, via the equivalent 2D cBD field theory. 
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The previous calculations in the cBD model were based on special analyticity conjectures. 
This is why we check our results for the asymmetric Coulomb gas, in the small-/? limit 
and close to the collapse f3 — 4 point. The small-/? series expansions are generated using 
the renormalized Mayer expansion in density, which is a simplified Coulomb version of the 
Feynman diagrammatic perturbation technique in 2D field theories. On the other hand, 
the important check of the results close to the opposite collapse point, based on screening 
properties of the Coulomb plasma, is original and has no counterpart in the field theory. 

The asymmetry in the strength of the particle charges brings two fundamental mod- 
ifications in the statistical behavior of the 2D Coulomb gas in comparison with its sym- 
metric version. Firstly, the large-distance exponential decays of the charge and density 
correlation functions are characterized by the same correlation length (in the symmetric 
Coulomb gas [19, 20], the charge correlation length is twice larger than the density one). 
Such behavior does not occur in the Debye-Huckel treatment of the model. Secondly, the 
mass of the lightest particle in the spectrum of the cBD model, which determines the 
thermodynamics and the asymptotics of the particle correlations, diverges at the collapse 
point (5 — 4. As a consequence, the truncated particle distributions are trivially equal 
to zero at (5 — 4. This is in contrast to the symmetric 2D Coulomb gas with finite and 
nonzero particle distributions at its collapse free-fermion point. 

The other cases of the asymmetric Coulomb gases do not belong to the family of 
integrable 2D Toda field theories. On the other hand, the 2D OCP, which corresponds to 
the extreme charge-asymmetric case, is exactly solvable at its free-fermion point [11, 12]. 
The 2D OCP has a field representation [48] which resembles the one of the 2D complex 
Liouville model with a kind of "background" charge. Although this theory does not 
belong to the Toda theories, it is integrable at least at the aforementioned free-fermion 
point. It might be therefore useful to explore its integrability properties at an arbitrary 
temperature, first in the classical limit and subsequently at the quantum level. 
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Appendix A: Small-/? expansions 

Let us consider a general fluid composed of distinct species of particles {a} with the cor- 
responding position-dependent densities {n(r, a)}. The particles i and j interact through 
the pair potential v(i,ai\j,<Tj), where vector position r-j is represented simply by i. The 
technique of the bond-renormalization in the ordinary Mayer expansion in density (for 
details, see Refs. [16, 31, 32]) is based on an expansion of each Mayer function in the 
inverse temperature /3, and on a consequent series resummation of two-coordinated field 
circles. The renormalized i^-bonds are given by 
K -0v 

o~v~~o = o- — -o + o- • -o + o- • • -o + ... 

l,<7i 2,cr 2 l,o"i 2,cr 2 l,en 2,cr 2 l,0i 2,a 2 

or, algebraically, 

ir(l,0i|2,0 2 ) = [-/3u(l,0i|2,0 2 )] 

+ E / d3 [-^(l,0i|3,0 3 )] n(3,<7 3 ) K(3,a 3 \2,a 2 ) (A.l) 

It is straightforward to verify by variation of (A.l) that it holds 

^ A (1 '" l|2 f ) = K(l, ^3, a 3 )K(3, a 3 |2, a 2 ) (A.2) 
dn(3,a 3 ) 

The excess Helmholtz free energy F ex is expressible in the renormalized format as follows 



-(3F cx [n\ = • • + D<°)[ n ] + J^D^n] (A.3a) 



8=1 



where 



d<® = <- ^> + / \ + ; ; + ... (A.3b) 

* * *- - - -4 



is the sum of all ring diagrams which cannot undertake the bond-renormalization proce- 
dure and 





etc. (A.3c) 



£)(1) £>(2) £)(3) 21 £(4) £)(5) 



are all remaining completely renormalized graphs. The free energy is the generating 
functional for one-body, two-body, etc. densities. The density-fugacity relationship is 
generated via 



In 



21(1,01) 



6(-(3F e *) 
571(1,01) 



o 



s=l 

where d^°'(l, a{) = 5D^ ' j Sn(l, a±) can be readily obtained in the form 

1 



(A.4a) 



en) = ^ lim [K(l, 0i|2, 2 ) + /^(l, 0-i|2, a 2 )] 



and 



5n(l, 0-1) 

The direct correlation function is generated via 

5 2 (-/3F ex ) 



(A.4b) 



(A.4c) 



c(l, 0"! 2,0-2 



5n(l, 0i)5n(2, 2 ) 

00 

0-----0 + c (o) (l,0i|2,0 2 ) + 5]c (s) (l,0i|2,0 2 ) (A.5a) 

1,01 2,0 2 



s=l 



where c^(l, 0i|2, 2 ) = <5 2 -D^°' ) /^(l, 0i)^(2, 2 ) corresponds to the Meeron graph, 

C (°)(l,0i|2,0 2 ) = 1 <3^ P a = ^ 2 (1, ^12,02) (A.5b) 

and 



c w (l,^i|2,02^ 



(A.5c) 



d>n(l, 0i)5n(2, a 2 ) 

The pair correlation function h, defined by relations (2.13) - (2.16), is related to c via the 
Ornstein-Zernike (OZ) equation 

/i(l,0i|2,0 2 ) = c(l,0i|2,0 2 ) 

+ ]T / d3 c(l,0i|3,0 3 ) n(3,0 3 ) /i(3,0 3 |2,0 2 ) (A.6) 

Notice that the functional derivatives with respect to the density field generate root circles 
not only at obvious field-circle positions, but also on K bonds according to formula (4.2), 
causing their right K-K division. 
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Let us return to the 2D asymmetric Coulomb gas with two kinds of particles a = 1, 2 
of charges (qi — l,q 2 — —1/2), interacting via the logarithmic interaction 



v(i,<Ti\j,<Tj) = qa t q* 3 v(i,j) 
v(i,j) = — In | z — j j 



(A. 7a) 
(A.7b) 



We consider the infinite-volume limit, characterized by homogeneous densities n(l,a) 
n a constrained by the neutrality condition J2a Q.c n a = 0, so that 



n 

ni= 3' 



2n 

n 2 = - 



(A.8) 



with n being the total particle density. Two-body functions are both isotropic and trans- 
lationally invariant, c(l, cr x j2, a 2 ) = c ai(72 (\l — 2|), /i(l, cr x j2, cr 2 ) = h ai<T2 (\l — 2|). It follows 
from Eq. (A.l) that the renormalized bonds exhibit the same charge- dependence as the 
interaction under consideration (A. 7), 



K{\a x \%a 2 ) = q ai q a2 K(l,2) 
K(l,2) = -(3K ( K \l-2\ 

Here, K is the modified Bessel function of second kind and 



K 



2irp m + 



n 2 



1/2 



\Jirf3n 



is the inverse Debye length. 



(A.9a) 
(A.9b) 



(A.10) 



A.l Density-fugacity relationship 

The small- x expansion of the modified Bessel function K (x) reads [44] 



oo x 2j 



oo x 2j 



K (x) = - In (-J g 2277JIP 4 ^ ^777m ; '(•>' - 1 ; 



(A.ll) 



o2*W 

where ^(l) = — C. Thus, in the uniform regime with the interaction (A. 7) and the 
renormalized interaction (A. 9), Eq. (A. 4b) yields 



C + In 



K 



a =1,2 



(A.12) 
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The first term on the rhs of (A. 4a) disappears due to the charge neutrality. Consequently, 



In 



C 



(7 = 1,2 



(A.13) 



where V — > oo is the volume of the system. For the dimensionless quantity of interest 
ra 1 -^ 8 /^) 1 / 3 , (A.13) gives 



1-/3/8 



(^f) 1 ^ 2 2 /3 



cxp 



2C+ ln[^ 



8 \ 3 <9ni ' 3 <9n 2 



19 2d 



E 

s=l 



(A. 14) 



When all {D^/V} are set to zero, we have the Debye-Hiickel approximation valid in 
the /? — * limit and, indeed, (A. 14) then reproduces the leading term of the exact (3- 
expansion (3.15). The scaling form of the renormalized interaction bonds (A. 9) permits us 
to perform the /^-classification of D^/V. Let the given completely renormalized diagram 
be composed of N s skeleton vertices and L s bonds. Every dressed bond K brings the 
factor —(3 and enforces the substitution r' = nr which manifests itself as the factor 1/k 2 
for each field-circle integration ~ J rdr. Since there are (N s — 1) independent field-circle 
integrations in we conclude that 

mm = ^.-^mizil (A . 15) 

In the sketch (A. 3c), the only diagram which contributes to the f3 2 order is the next 
diagrams D^ 2 \ and constitute the complete set of contributions to the (3 3 

order, etc. 

We are interested in the lowest correction to the Debye-Hiickel limit, so only the 
diagram has to be analyzed. The contribution of this diagram is expressible as 
D^(n u n 2 ) 



V 



^(^-2m| + |)/d 2 r[-^ Mr 



1 _^(ni-n 2 /8) 2 fd 2 r 



-KfJr) 



(A.16) 



2\3V ni+n 2 /4 J 2n 
To evaluate the integral in (A.16), we can make use of the 2D Fourier components of 
K (r) and Kq{t) (see Appendix of Ref. [16]) to derive 

In [(k/2) + + (k/2f 



J 27T ° V J (l + k 2 )Jl 



+ {k/2f 



(A.17) 
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After the substitution k = 2sin(t/2), one gets 
/v,-;(r) = - 

1 3=1 



2tt 



1 ^ sin(7rj/3) 
v^fe J 2 



1 

72 



*(§)+*(§) 



- * 1 1) - * (I 



where ^//(rr) is given by (3.16). Using the readily derivable relationships [44] 

'x + r 



i>'{x) = 



V/ 



X 



lfj'( X ) = ip'(x + 1) + 



1 



X/ 



valid for any x, one finds 



d 2 r 



1 
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3^' 



2^ 



3^' ( § ) " 2tt 2 



/? 2 



1728 



(A.18) 



(A.19a) 
(A.19b) 
(A. 19c) 



(A.20) 



(A.21) 



The consideration of (A.20) in (A. 16) leads to 

(id 2 a \ d<v 

\ 3 drii 3 dn 2 ) V 

where (A. 8) was taken into account. Finally, inserting (A.21) into (A. 14), the term of 
the order f3 2 is reproduced exactly in the exponential of the small-/? expansion (3.15), 
generated from the exact density-fugacity relationship (3.14). 



A. 2 Large-distance asymptotics of particle correlations 

In the Fourier picture, the OZ equation (A. 6) reads 



h a a'(k) = c aa >(k) + 2ir ^ c aa »(k)n a »h a » a >(k) 

a" =1,2 



(A.22) 



It follows from (A. 5) that, at the lowest order in (3 (Debye-Huckel approximation), 



w(r) 

Caa'(k) 



f3q a q al lnr 
1 



k 2 



(A.23a) 
(A.23b) 
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Consequently, from (A. 22), 



h aa >(k) 
h a a'{r) 

and the asymptotic form of h aa '(r) is 



k 2 + k 2 
-Pq a q a >K (Kr) 



(A.24a) 
(A.24b) 



/ 7T \V 2 

h a a'(r) ~ -fiqaQa' ( ^ J exp(-«r) 



(A.25) 



reproducing thus the expected result (4.9) with A CTO -/ = (3q a q a > + 0(/3 2 ) (4.13a) and m = 
k[1 + 0(/3)] (4.13b). 

With the aid of Eq. (A. 5), the direct correlation function is determined up to the f3 2 
order as follows 



Caa'(r) = Pqaq*> lnr + —qlq^K^Kr) 



(A.26) 



Taking k 1 = l/y/ix(3n as the unit of length, one has in the Fourier space 

\k/2) + yf\ + (fc/2) 



1 /3 2 



In 



fc Jl + (A;/2) 2 



(A.27) 



The OZ equation (A. 22) implies 



hn(k) = 


ii 

d 1 


f „ /3 2 (3 + 4A; 2 ) 


h 12 (k) = 


ii 

d 1 


f* ln[ 

[2 4^4 + fc 2 [ 


h 2 2(k) = 


ii 

d I 


f /? ! /3 2 (6 + A; 2 ) 
[ 4 16/cV4 + k 2 



In 



A; + V / 4TI 2 



In 



(A.28a) 
(A.28b) 
(A.28c) 



where the denominator is 

d = 1 + k 2 - 



In 



k + + W 



(A.29) 



4kV4 + ¥ 

The asymptotic behavior of h aa >(r) is governed by the poles of h a(7 >(k) closest to the real 
axis. When (3 — > 0, there poles are at k = ±i (or k 2 = — 1). In close analogy with Ref. 
[19], to find the /3-correction of these poles, it is sufficient to expand the /^-dependent part 
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of the denominator d around say k = i up to the first order in (k — i): 

7T(3 \(3 



d = l + k 2 + 



— (9 + 8v / 37r)(A;-i) + --- 



= (k 2 + m 2 ) 
= (k 2 + m 2 ) 



24x/3 108 

i(3 (9 + 8 v / 37r)(A;-i) 



108 



k 2 + m 2 



+ 



where 



m — 1 + 



48^ 



The prefactors X aa >, defined by 



Ka'{k) = -X 



+ 0(J3*) 



1 

k 2 + m 2 



are deducible from (A. 28) (with numerators evaluated at k 2 = — 1) and (A. 30), 



Al2 — 

A 22 = v + '— \ ^7 



24 72^/ 
2 2 1 24 36^/ 



0(/3 3 ) 



(A.30) 
(A.31) 

(A.32) 

(A.33a) 
(A.33b) 
(A.33c) 



4 4 \24 72^/ 

The large-distance behavior of /w( r ) is of type (4.9): the prefactors X aa > (A. 33) and 
the parameter m (A.31) have the small-/? expansions (4.13a) and (4.13b), respectively, 
confirming thus the predictions of the form-factor method. 



Appendix B: Collapse point 
B.l Density- fugacity relationship 

Based on the analysis in Section 2 applied to our specific Q = 2 case, close to the collapse 
point j3 — 4 (b — 1), the two-body densities behave like 

nia(r) ~ z x z 2 ^l 2 )\ h ^l 2 (B.la) 
nu(r) oc r 1 (B.lb) 
n 22 (r) oc r _1 (B.lc) 
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in the short-distance limit r — > 0. The average in (B.la) is taken over the cBD action 
(3.4) with b -> 1. 

The one-body and two-body densities in a general Coulomb system satisfy the elec- 
troneutrality sum rule [33] 

-.A^v/d'-vW (B.2) 

where a numerates the charged species. For our asymmetric 1/ — \ 2D Coulomb gas, one 
has in particular 

~ = J d 2 r nu(r) d 2 r n i2 (r) (B.3a) 

| = |d 2 rn 12 (r)-i|d 2 rn 22 (r) (B.3b) 

It follows from (3.17) that, for fixed zi and z 2 , the density n exhibits the singularity of 
type (1 — (3/A)~ 2 ^ 3 as (3 — > 4 _ . This singularity originates in (B.3) as a result of the 
short-distance integration over the two-body densities ni 2 (r) oc r _/3 / 2 ; neither nu oc r 1 
nor n 22 (r) oc r _1 give a diverging contribution as (3 — > 4~ after being integrated out over 
short distances r. The problem with Eqs. (B.3) is that they imply the dependences of n 
on the relevant (diverging) integral / d 2 rn i2 (r) with two different prefactors. We have to 
admit that the collapse mechanism alters the form of the sum rules (B.3). The particles 
can be divided into three basic groups: one third of particles is of type 1 (charge q\ — +1 
and density n\ = 1/3), one third of particles is of type 2 (charge q 2 = —1/2 and density 
n 2 = 1/3) and the remaining third of particles of type 2' (charge qy = —1/2 and density 
n 2 ' = 1/3) is excluded from the collapse phenomenon (they only feel the charge +1/2 of 
each collapsed pair of 1,2-particles and do not enter into diverging screening integrals). 
Under this assumption, close to (3 — 4, the electroneutrality sum rules (B.2) are modified 
as follows 

-q 2 >n 2 > - q a n a ~ q a * j d 2 r n aa *(r), (7 = 1,2 (B.4) 
where 1* = 2 and 2* = 1. The couple of Eqs. (B.4) is consistent and implies the only 
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singular relationship 

n ~ 3 J d 2 r n 12 {r) as (3 -> 4" (B.5) 

We now put the short- distance expansion of ni2(V), Eq. (B.la), into the integral on the 
rhs of (B.5) cut at some finite |r| = L (L is a length over which the Coulomb interaction 
is screened) and obtain 

3-7T 

n ~ z 1 z 2 (^ 2 )\^ 1 — — (B.6) 

To evaluate the mean value of the exponential field, we apply the conjectured formula 
(3.8) for a = 1/2, b — > 1 and £ = 1. Consequently, 



<e^ 2 >| 



b-»i 



x/2 



* 2 r(3/4) / 
*r(i/4)l 4, 



1/3 



dt 

X< ' XP Wo 27 



-2t 



fB 7") 

_cosht(2cosht - 1) ~ J J 1 ' ; 

where we have applied F(x + 1) = xT(x). It is simple to show by using the integral 
representation of the logarithm of the Gamma function (3.11) that 



dt 

rX1, Wo 2t 



1 



-2t 



cosht(2 cosht — 1) 



<2\V 2 T(3/4) 
(2tt) 2 V3/ f(l74) 



(B.8) 



Eqs. (B.6) - (B.8) reproduce exactly the collapse singularity (3.17) deduced from the 
exact density-fugacity relationship (3.14). 



B.2 Thermodynamics 

Although the thermodynamics of the system close to the collapse point is derivable directly 
from the collapse n — z relationship discussed in the previous subsection, there exists 
another simpler derivation in the spirit of an independent-pair approximation by Hauge 
and Hemmer [6]. As was already mentioned, close to the collapse point, from the total 
number of iV particles, iV/3 particles of type 1 and iV/3 particles of type 2 form pairs. 
Their statistical weights contribute dominantly to the configuration integral 







N/3 r 
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/ L d 2 r r-^ 2 
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.£2(1-/9/4) 



N/3 



as (5 -> A' 



;b.9) 



Thence 

S = ^W = iFW I(w7i) +0(1) (R10) 

in full agreement with the Laurent series (3.24). 
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